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ABSTRACT 

Right handed neutrinos with mass of O (10 12 — 10 13 ) GeV are required to imple- 
ment the see-saw mechanism and generate neutrino masses capable of playing a role 
in structure formation. Moreover models of fermion masses often relate the Yukawa 
couplings involving these neutrinos to the up-quark Yukawa couplings. Here we study 
the effects of such couplings on the radiative corrections to quark masses. We find that 
b — t equality at Mqijt may still give the correct mi/m T — ratio at low energies, but 
only if there is large /i — r mixing in the charged leptonic sector. We propose specific 
mass matrix "textures" dictated by a U(l) family symmetry whose structure preserves 
nib = n^r at Mqut- In these schemes, due to the large — v T mixing, it is possible to 
give a simultaneous solution to the solar neutrino deficit and the atmospheric neutrino 
problem. 



* SERC Senior Fellow. 



1 Introduction 



Although the Standard Model is remarkably successful in describing strong and elec- 
troweak phenomena it still leaves many questions unanswered. In particular one would 
like to understand the origin of fermion masses and mixing angles which in the Standard 
Model appear as arbitrary parameters. An obvious possibility is that some structure 
additional to that of the Standard Model is responsible for the pattern of masses and 
mixings that we see at low energies. Support for such a stage of unification has been 
obtained for, although unification |ffl, ^] on its own does not agree with experiment, 
when combined with supersymmetry it leads to very successful predictions || for the 
gauge couplings, the pattern and magnitude of spontaneous symmetry breaking at the 
elecroweak scale [Q] and the bottom - tau (b - r) unification |5|. A further indication 
that additional symmetries beyond the Standard Model exist, has been the observation 
that the fermion mixing angles and masses have values consistent with the appearance 
of "texture" zeros in the mass matrix || [|, ||, |9|, [H], 1 1 , [l2j . Such zeros may indicate the 



presence of an additional broken symmetry. When unbroken only the third generation 
is massive and all mixing angles are zero. However, symmetry breaking terms gradually 
fill in the mass matrix and generate a hierarchy of mass scales and mixing angles. 

In this paper we will consider the implications for the neutrino sector extending the 
analysis presented in [11]. We shall consider models with right-handed neutrinos in 



which both Dirac and Majorana masses for the neutrinos are present, and the "see-saw" 
mechanism automatically explains the lightness of neutrinos relative to the charged 
fermions. 

Neutrino data from various experiments can be explained if there is mixing between 
the various types of neutrinos. The solar neutrino puzzle can be resolved through 
matter enhanced oscillations (preferably between v e and states) with a mixing angle 
somewhat smaller than the 1/3 of the corresponding Cabbibo angle of the quark sector, 
Vq Km . For this explanation to work, the squared mass difference of the two types of 
neutrinos involved in this phenomenon should lie in a very narrow region. The specific 
ranges for the angle and mass squared given by the latest experimental data are: 

For matter enhanced oscillations in the sun: 

sin20 ex = (0.39 - 1.6) x 10~ 2 , Am 2 = (0.6 - 1.2) x lO^eF 2 , (1) 

and for vacuum oscillations: 

sm28 ex > 0.75, Am 2 = (0.5 - 1.1) x 10~ 5 eV 2 . (2) 

If we wish to avoid fine tuning problems, it seems necessary to assume that such small 
differences in neutrino masses can be obtained only if the v e ^ neutrino masses them- 
selves are also of the same order. Finally, if neutrinos play a role in structure formation, 
providing a hot dark matter component, then the heavier neutrino(s) should have mass 
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in the range ~ (1 — 6) eV, the precise value depending on the number of neutrinos that 
have masses of this order of magnitude. 

What symmetry could explain this pattern of masses? If gauge symmetries have some- 
thing to do with the hierarchical fermion mass spectrum a similar hierarchy may be 
expected to hold for the unknown neutrino masses too. In a previous study it was found 
that the observed hierarchical mass spectrum of the charged fermions (quarks and lep- 
tons) follows naturally if we extend the gauge group of the minimal supersymmetric 
standard model by a Ufd{^-) family type symmetry. The extension of this model to 
include the right handed neutrino in the theory resulted in a similar structure of the 
neutrino sector as well leading to the following general structures [fj"l]1 : 

1. ) The solar neutrino puzzle can be explained via v e — ► i/„ oscillations. The hierarchi- 
cal mass spectrum leads to the conclusion that m Uft ~ \J^ m j l ,e while m Ue <C m Vji . This 
also fixes the right handed neutrino scale Mjv through the effective Major ana mass 
matrix resulting from the usual "see - saw" mechanism 

< ff = -\^v D M N l m T VD (3) 

2. ) The simultaneous solution of the solar neutrino problem and the interpretation 
of the v T mass as a hot dark matter component through the effective light Majorana 
neutrino mass matrix, requires a right handed neutrino scale of the order of M VR ~ 
10 12 - 10 13 GeV. 

3. ) There is no natural solution of the atmospheric neutrino problem unless a consid- 
erable fine-tuning of the coefficients in the neutrino mass textures occurs. This follows 
because the U(l) symmetry that was used to derive the above textures together with 
simple spontaneous breaking gives only a hierarchical mass spectrum and small mixing 
angles for all fermion mass matrices. 

One additional implication of the structure emerging from the U{1)fd symmetry is 
that the right handed neutrinos have Yukawa couplings of the same order as the up 
quarks. This in turn affects the radiative corrections in the model and in particular 
the expectations for gauge unification and for the mj/m T ratio. The implications of 
such large couplings has already been explored in refs.p3|. Here we develop these 
analyses in two respects. Firstly we present a semi-analytic analysis of the radiative 
corrections that allows us to analyse the possibilities for maintaining b — r equality at 
the GUT scale even in the presence of these radiative corrections through large fj, — r 
mixing giving a mechanism to evade conclusion 3.) aboveQ. This leads us to consider 
schemes based on the U(l) family symmetry which naturally generate such mixing. 
The implications for the neutrino mass spectrum are then explored. In addition the 
semi-analytic approach is supported by a full numerical calculation. In section 2 we 

^or another mechanism evading this result through the introduction of spontaneous U(l) 
breaking via several Higgs scalars see ref . Q . 
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give a semi-analytical approach to the renormalisation group equations, in the presence 
of right handed neutrinos. These equations are used in section 3 in order to get some 
direct intuition on the effects of the heavy neutrinos. The explicit form of the solutions 
makes it easy to see how b — r equality at a GUT scale may be made consistent with 
the parameter spectrum at low energies, by sufficient \x — r mixing in the charged 
leptonic sector and for a relatively heavy strange quark. This is discussed in section 
4. In section 5 we give the resulting predictions for the heavy and light Majorana 
neutrino mass matrices and eigenvalues, with a mixing which is of the correct order 
of magnitude in order to explain the atmospheric neutrino problem. In section 6 we 
present a numerical approach to the renormalisation group equations, which depicts 
the observations of section 2. Finally, in section 7 we give a summary of our results. 



2 RGE with RH-neutrinos: a semi — analytic 
approach 

From the above it is clear that the interpretation of many important experimental facts 
is based on the existence of the right - handed partners vr. of the three left - handed 
neutrinos, where the scale of mass of these particles is at least three orders of magnitude 
smaller than the gauge unification scale, Mjj. Thus the running from the Unification 
scale, Mjj ~ 10 16 GeV, down to the scale of M UR , must include radiative corrections 
from vr neutrinos. After that scale, vrs decouple from the spectrum, and an effective 
see - saw mechanism is operative, c.f. eq( |3|). 

In the presence of the right handed neutrino, the renormalization group equations for 
the Yukawa couplings at the one-loop level are 

167r2 Jt ku = ( 3huh u + hoh D + 1 ■ Tr[h N h) N ] + / • Tr[3h v h\j] - I ■ G v ) h v , (4) 

16n2 dt hD = ( 3hok D + h u h u + 1 ■ Tr[3h D h j D ] + I ■ Tr[h E h ] E ] - I ■ G D ) h D , (5) 

167r2 Jt kE = ( 3hEk E + h * h N + 1 ■ Tr l h EhU + I ■ Tr[3h D h j D ] -I-G E ) h E , (6) 

l ^ 2 Jt kN = ( h E h E + 3h Nh^ + I-Tr[3h u h\ I ) + I-Tr[hW N }-I-G N ' S )h N .(7) 

where h a , a = U, D, E, N , represent the 3 ® 3 Yukawa matrices for the up and down 
quarks, charged lepton and Dirac neutrinos, while I is the 3 <8> 3 identity matrix. Fi- 
nally, G a = J2i=i c a9i(t) 2 are functions which depend on the gauge couplings with the 



coefficients c^'s given by [15, 13]. 



{4/};=i, 2 ,3 = {^' 3 'y}' { c £>}i=l,2,3 = Irg' 3 ' y}> ( 8 ) 
{<4}i=iA3 = {|,3,oj , {4} i= i l2 ,3 = {|3,0}. (9) 



3 



Consider initially the simple case where only the top and Dirac - type neutrino Yukawa 
couplings are large at the GUT scale (i.e. the case of small tanj3 scenario). Let us start 
assuming that the top and neutrino Yukawa couplings are equal at the Unification scale, 
ht(Mu) = his[{Mu)i a relation which arises naturally not only in our case but in most 
of the Grand Unified Models which predict the existence of the right handed neutrino. 



As in the case of the charged fermions, we will consider only hierarchical textures [11] 
for the right handed neutrino Majorana mass matrices, i.e. M Vl <C M U2 <C M V3 . If we 
work in a diagonal basis we can considerably simplify the above equations. Then, for 
the range Mjj to , the renormalization group evolution of the Yukawa couplings of 
third generation, can be written as follows 

^ 2 jh t = (6h 2 t + h% - Gu) h t , (10) 
16tt 2 -/^v = ^h 2 N + 3h 2 -G N )h N , (11) 

167r2 Jt h = ( h t- G D) h b, (12) 
l^ 2 j t K = (h 2 N -G E )h T , (13) 

Below Mjy, the right handed neutrino decouples from the massless spectrum and we 
are left with the standard spectrum of the MSSM. For scales Q < Mjy the gauge and 
Yukawa couplings evolve according to the standard renormalisation group equations. 
In addition, the effective Yukawa coupling of the light neutrino mass matrix (Q) evolves 



according to [16| 



167r 2 -^ =h v (I-Tr[6h u hl J }-G u ) + h v h E h j E + h j E h E h u (14) 

with G v = 2g\ + §g\. In order to gain an insight into the effects of new couplings 
associated with the vr in the renormalisation group running we integrate the above 
equations in the region < Q < Mjj. We denote the top and ur Yukawas by ha 
at the unification scale, while bottom and r are denoted with hb , h TQ respectively. We 
get 

h t (t) = Tu(t)hG$tN (15) 
M*) = 7N(t)h G £& (16) 
h(t) = 7d(0K6 (17) 
h T {t) = 7E(t)h T0 £ N (18) 

where the functions J a (t) depend purely on gauge coupling constants and are given by 
7a (t) = exp(^-2 J t G a (t)dt) (19) 

- ii (r) 



4 



= ^..jg^zMj^, (21) 

The £j's (i = t, iV, 6, r) are given by the integrals 

Si = ^Ttbil^ (22) 

The values of the parameters & can be determined at any scale by numerically solving 
the renormalization group equations. As a general remark, we note that the initial 
condition for £j is ii{tjj) = 1, while at any lower scale Q < My, £i(Q) < 1. 



3 Heavy neutrino effects : an insight 

We start by investigating the b — r Yukawa coupling solutions. Thus, in the case of 
small tan(3, we can relate their solutions at the scale Mn in terms of the initial values, 
by the following equation 

h b (t N ) = p^ t —h T (t N ) (23) 

IE 

with p = h ^ b £ N ■ In the case of b — r unification at My , we have h T0 = hb , while in 
the absence of the right - handed neutrino £jv = 1, thus p = 1 and the m& mass has 
the phenomenologically reasonable prediction at low energies given by the approximate 
one-loop formula 

m b = r) b £t — m T (24) 

IE 

where % is the renormalization group coefficient in the m t — range. In the presence 
of vr however, if h To = hj, at the GUT scale, the parameter p is no longer equal to 
unity since £jv < 1. In fact the parameter £jv becomes smaller for lower M/y scales. 
Therefore in order to restore the correct mb/m T prediction at low energies we need 
p = 1 corresponding to 

hb = KAn (25) 
Hence it is obvious that we need a r— Yukawa coupling h T0 , larger than hb at My to 
compensate for the factor £jv arising from the presence of vr. 

It is interesting that £/v can be given in this case in terms of the values of gauge and 
Yukawa ratios at only, irrespective of the initial conditions 

t hb o 

= h b j E , h£YN yi/3 / 2g x 

On the other hand, the top mass at the scale = In(M^) can also be expressed 
formally as follows 

m t (t N ) = h G iu(t N )$(t N )£ N (t N )-j=sin(3 (27) 
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where v = 246GeV and tanj3 is the ratio of the two Higgs vev's. Again, in the absence 
of ur this reduces to the well known one loop approximate formula which coincides 
with the above for £jy = 1. In the present case however, this prediction corresponds 
effectively to a smaller initial value of the top Yukawa coupling of the order 

h' G = h G ^N(t N ) (28) 
For Kg > 1) however, due to the infrared fixed point property of the top - Yukawa 



solution [17], the nit - prediction is not going to alter significantly. For the same 
tan/3, one will get almost the previous top mass prediction, reduced at most by 2%. In 
contrast, in the small tan(3 scenario where hb *C 1, one naturally expects that bottom 
mass will be affected by the presence of vr. For ~ 10 13 GeU for example and 
he > 1; we can estimate that £(iA?) ~ 0.89 thus, there is a corresponding ~ 10% 
deviation of the r — b equality at the GUT scale. 

Furthermore, we have seen that in order to recover the correct m^/m-,- relation at low 
energies, it is necessary that hf,fi/h Tt o < 1 as long as M/v < M\j. This can be done 
in two ways: Either we can keep the same value of the b - Yukawa and increase the 
r- Yukawa by a factor ^j- 1 , or decrease the bottom coupling by a factor £jy in the first 
case, the angle f3 remains the same and the top mass unaffected. In the second case, 
in order to retain the same absolute value of the initial parameters for the 6, r masses 
we need to increase cos(5. This results to a corresponding decrease of the top mass 
prediction. 

We will present a detailed numerical analysis of the above in section 5, where two loop 
effects from the gauge couplings are taken into account. In the next section we first 
propose a scheme in which the bottom-tau unification is restored. We will show that 
this may result in a solution of the solar neutrino deficit and the atmospheric neutrino 
problem. 



4 Restoration of bottom — tau unification 

Given the results of section 3, it is natural to ask if Grand Unified models which predict 
the b — t equality at the Unification scale, exclude the experimentally required and 
cosmologically interesting region for the neutrino masses. To answer this question, we 
should first recall that the b — r - equality at the GUT scale refers to the (33) entries 
of the corresponding charged lepton and down quark mass matrices. The detailed 
structure of the mass matrices is not predicted, at least by the Grand Unified Group 
itself, unless additional structure is imposed. It is possible then to assume (m^)^ = 
( m Z))33 an d a specific structure of the corresponding mass matrices such that after the 
diagonalisation at the GUT scale, the {in 5 E )^ and (771^)33 entries are no-longer equal. 

To illustrate this point, let us present here a simple 2x2 example. Assume a diagonal 
form of mP D at the GUT scale , m° D = diagonal(cmo,mo), while the corresponding 
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entries of charged lepton mass matrix have the form 




m (29) 

These forms of m Q D ,m Q E ensure that at the GUT scale (m^-,)^ = (7740)33. However, at 
the low energies one should diagonalize the renormalised Yukawa matrices to obtain the 
correct eigenmasses. Equivalently, one can diagonalise the quark and charged lepton 
Yukawa matrices at the GUT scale and evolve separately the eigenstates and the mixing 
angles. Since m a D has been chosen diagonal there is no need of diagonalization and the 
mass eigenstates which are to be identified with the s,b - quark masses at low energies 
are given by 

m s = cy D m , m b = 7rjm 6 (30) 

with mo = hb -^cos(3. To find the charged lepton mass eigenstates we need first to 
diagonalise m E at Mqut- We can obtain the following relations between the entries 
e, d of mg and the mass eigenstates m^,m" at the GUT scale. 

m2. — m?. 

d=(— £-1) (31) 

m 

n m?. m° 

e 2= + l F^-l (32) 

In the presence of right handed neutrino, the evolution of the above r— eigenstate 
down to low energies is that described by (|l^) with m T0 = h TQ -^cos(3. By simple 
comparison of the obtained formulae, we conclude that, to obtain the correct m T /mi, 
ratio at mw while preserving the b — r unification at mcuT, the m E entries should 
satisfy the following relations 

' J2 




- 1) = e 2 (33) 

The above result deserves some discussion. Firstly we see that it is possible to preserve 
b — t unification by assuming 2 — 3 generation mixing in the lepton sector, even if the 
effects of the vr states are included. Secondly, this mixing is related to a very simple 
parameter which depends only on the scale Mn and the initial condition. The range 
of the coefficient c in the diagonal form of the mP D - matrix, can also be estimated using 
the experimental values of the quark masses m s ,mb. An interesting observation is that 
the usual GUT - relation for the (22) - matrix elements of the charged lepton and 
down quark mass matrices, i.e., {rriE)<i2 = — 3(mo)22, which in our case is satisfied for 
c = —d/3, implies here a relatively heavy strange quark mass m s ~ 200MeV. Smaller 
m s values are obtained if — 3c/ d < 1. 



2 An alternative mechanism which restores the correct mb/m T ratio in the presence of vr 
was proposed in Q. 
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We turn now to a consideration whether the hierarchical structure of the lepton mass 



matrix corresponding to eq (|33|) can be obtained by a simple U(l) symmetry |10 |, 



In |10[ a viable fermion mass matrix structure was constructed following from a spon- 
taneously broken U(l) gauge symmetry. In this the form of the down mass matrix 
is 



(We have concentrated for simplicity in the case of the two heavier generations which 
are relevant here.) Note that we have suppressed all Yukawa couplings in writing this 
mass matrix - only the order of the matrix elements allowed by the broken symmetry 
are given. These Yukawa couplings are assumed to be of order 1 and the object of the 
exercise is to demonstrate that the hierarchical structure of the fermion masses may 
come only from symmetry considerations. Eq(|34]) is diagonalised by the orthogonal 
matrix 

H-'-O (35) 

where e = 0.23, in order to fit the down quark masses and mixing angles. 

The structure of the lepton mass matrix following from the U(l) symmetry (again for 
the heavier generations) is 

W 1 ) ^ 

where = 1-6= f^E^-. and in @, Q the cases (3 = 1/2 and = 1 had been 




considered. Both possibilities are in very good agreement with the measured masses and 
mixing angles. The important fact here is that j3 can be determined by the requirement 
that b — t mass ratio be correctly given when heavy neutrinos, which become massive 
at an intermediate scale M/v < My, are present. Allowing for the unknown coefficients 
of 0(1) we see (cf. eqs([3rj|) and (^)) that both j3 = 1/2 and (5 = 1 are in reasonable 
agreement with the above expectation^]. Now the diagonalising matrix is given by 



Vl-e 

Obviously, there is a large mixing in the 2 — 3 lepton sector of the obtained solution 
which may lead to interesting effects in the rare processes like r — > /U7 and neutrino 
oscillations. In the simplest realisation of this scheme we expect hb ~ ht because in the 
limit e, e — > the U(l) quantum numbers of the light Higgs Hi 2 allow them to couple 
to the third generation and a SU(2)i®SU(2)ji symmetry of the couplings ensures equal 
Yukawa coupling. Thus this model applies only to the large tan f3 regime. However if 



3 Here we assume the field spontaneously breaking U(l) carries half - integral U(l) charge 
so we do not have the Z 2 symmetry of |LQ] . 
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there is an additional heavy state, Hi, Hi, i = 1 or 2, with the same U(l) quantum 
number then mixing effects can generate different hf, and ht couplings allowing for any 
value of tan /3. 



5 The Effective Light Majorana Mass Matrix 

We have seen in section 3, that we can obtain with a U{\) family symmetry a charged 
lepton mass matrix with the required large mixing in the two heavier generations by 
choosing the one free parameter, b. The choice 6 = 1/2 gives a very good agreement with 
the charged lepton masses and the bottom-tau relation in the presence of vr with mass 
Mjv ~ 10 13 GeV. Our next step is to determine the Dirac and heavy Majorana mass 
matrices. The general form of the Dirac neutrino mass matrix for arbitrary a, [3 = 1 — b 



is given by [11] 



£ 2|3a+/3| e 3|a| £ |3a+/3| 
M D = | £ 3H e 2|/3| € \f)\ J (3 g) 
e \3a+f3\ e \f3\ 1 

The Majorana masses for the right - handed neutrinos are generated by terms of the 
form urvrYi, where £ is a singlet scalar field -invariant under the SU(3) <8> SU(2)l <8> 
U(l)y gauge group- with charge under the U(1)fd family symmetry. For the 
various choices of Qs we may then form the possible "textures" for the heavy Majorana 
mass matrix. For example, when Qy, = — 2a±, that is when the field £ has the same 
charge with the Higgs fields, only the (3, 3) entry of the mass matrix appears for an 
exact U(l) symmetry and is of order unity. 

However, at a later stage the Abelian symmetry is broken by standard model singlet 
fields 9 and 9 with U{\) charge ±1/2, which acquire vacuum expectation values along 
a D - flat direction. At this stage, invariant combinations involving the 9, 9 fields are 
generated, filling the rest of the entries in the mass matrices in terms of an expansion 
parameter fT0|| . Depending on which elements of the Majorana mass matrix we require 
to appear before the U(l) symmetry breaking, we can make six different choices of the 
charge Qy, which result to the M™ -?' - "textures" shown in Table 1. 

For a = 1, (3 = 1/2, we can obtain the specific forms for Dirac and Majorana textures 
compatible with the correct fermion mass predictions in the presence of the intermediate 
neutrino scale. In Table 2 we present the eigenvalues of the heavy Majorana mass matrix 
for this choice of a and j3. 

The analysis of the resulting mff* follows the same steps as in ref . [ 1 1 1 . In the matrices 

4 The full anomaly free Abelian group U(l) involves an additional family independent com- 
ponent U{1)fd- 
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/ g2|3a+/3| E 3\a\ ^3a+p\ ^ 

g|3a| E 2\P\ g|/3| 
\^ e|3a+/3| x J 






/ g3|2a+/3| g|3o+/3| g|3a+2/3| ^ 

g|3a+/3| g|/3| x 
y g|3a+2/3| x g|/3| j 






/ g2|3a+2/3| g|3«+2/3| g3|a+/3| \ 

g|3a+2/3| j 
K g3|a+/3| g|/3| g2|/3| J 






/ g|3a+/3| g|/3| x \ 

g|/3| g3|a+/3| g3|a+2/3| 
^ ! g3|a+2/3| g|3a+/3| y 






/ 1 g|3a+2/3| ^|3a+/3| \ 
g|3a+2/3| g2|3a+2/3| g3|2a+/3| 
g|3a+/3| E 3\2a+(3\ g2|3a+/3| 


) 




/ E \3a+2/3\ 1 \ 

! g|3a+2/3| g|3a+/3| 
^ gl/3| g|3«+/3| g|3a| y 





Table 1: General forms of heavy Majorana mass matrix textures. The specific 
textures of the text arise foro; = l, / 9 = l/2. 





( el ° \ 

e 2 

V i J 






/ e 15 \ 

-1 + e 
V 1 + e / 






r-j 






/ e 9 \ 
-1-e 2 

V i + e 2 y 






/ e 16 N 
e 14 

1 / 


1 




/ e 6 \ 

-1-e 2 
V 1 + e 2 y 





Table 2: Eigenvalues of Heavy Majorana mass matrix textures, for a — 1 and 
P = 1/2 





1 \ 

V i / 






/ e 25 \ 

e 9 

V 1/e / 




/ 
\ 


g 24 

e 8 


) 


/ 
\ 


/ e 33 

e 13 

1/e 7 , 


) 


/ e 40 ^ 
1/e 8 

V i/e 14 J 






/ e 32 \ 

e 6 

V 1/e 6 J 





Table 3: Eigenvalues of light Majorana mass matrix textures, for a — 1 and 
f3 = 1/2 
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we use 

e = e 1/2 , e = 0.23 (39) 

The eigenvalues of m e ff are given in Table 3. The order of the matrices in Tables 2 
and 3 corresponds to the one of Table 1. 

Note the interesting feature that the large mixing in the (2,3) entries of the charged 
leptons which we introduced in order to restore the b — r unification leads to a similar 
large mixing in the neutrino sector^]. This mixing is of the correct order of magnitude 
for a possible solution to the atmospheric neutrino problem. Indeed, the atmospheric 
neutrino problem may be explained in the case that large mixing and small mass split- 



ting involving the muon neutrino exists [IS]. Taking into account the bounds from 
accelerator and reactor disappearance experiments one finds that for v e — v„ or v T — v„ 
oscillations 

bml^ > 10- 2 eV 2 (40) 
sin 2 26^ a > 0.51 - 0.6 (41) 



In [11], such a large mixing was not present due to a residual discrete symmetry assumed 
for the 6 = 1/2 case. In that case, from the resulting ray* mass matrix we had been 
able to fit the COBE results and solve the solar neutrino problem (solution A). 

In the case of the large mixing discussed here we may also have a similtaneous solution 
to the solar and the atmospheric neutrino problems (solution B). However, the small 
mass splitting which is required between the neutrinos that mix in both the solar 
and atmospheric neutrino oscillations, make it impossible to account for the COBE 
measurements at the same time. This is a result of working in the minimal scheme 
with only one £ field present, which naturally leads to a large splitting between the 
neutrino masses. However, in the case that we add more singlets £ in the theory, it is 
possible to obtain a heavy Majorana mass matrix that leads to a solution of all three 
problems similtaneously |L"fl| . 

Going back to the case of a single £ field, whether we obtain the solution (A) or (B) 
depends on the predicted mass splitting between the two heavy neutrinos in each of 
the six choices of the heavy Majorana mass matrix. For a v T ~ 5 eV and x% = e 6 , e 8 
and e 10 , for i = 1,2,3, we obtain a muon neutrino mass m v = m UT Xi = 0.06,0.014 
and 0.003 eV respectively. This indicates that our solutions with a total splitting 
e 10 naturally lead to a solution of the COBE measurements and the solar neutrino 
problem. On the other hand, for m Ur rj 0.1 eV and x\ = e 6 , m u ^ = m Vr x\ = 0.0012 
eV, which may be marginally consistent with a solution to the atmospheric and solar 
neutrino problems (remember that coefficients of order unity have not yet been defined 
in the solutions). Since there are alternative schemes which lead to an explanation of 



5 The mixing in the (1,2) is of the C((^) 1 / 2 ) and negligible in (1,3). 
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the COBE measurements, other than hot and cold dark matter ^ we believe that the 
scheme (B) should be considered on equivalent grounds with the scheme (A). 



6 Numerical analysis 



In this section, we present the results of a numerical analysis of the effects of the 
heavy neutrinos in the renormalisation group analysis of masses, concentrating on its 
implications for lepton mass matrices with a large [i — r mixing. We start by giving a 
brief description of the procedure we are following. We compute numerically the low 
energy values of gauge and Yukawa couplings, taking into account two - loop effects of 
the gauge couplings, one loop RGEs for the Yukawas assuming an effective SUSY scale 
to account for low energy threshold effects. 

First, we check the procedure by reproducing the standard results when no right handed 
neutrino is present in the theory. We start at the unification scale Mjj using as inputs 
Mjj itself, the values of the common gauge coupling au, the top coupling ht G and 
hb 1 h TQ . In obtaining the low energy values of a ern , 03, and sin 2 6y/, we use the following 
ranges 

a em - 1 = 127.9 ± .1, a 3 = .12 ± .01, sin 2 w = .2319 ± 0.0004 (42) 
The top mass mt is obtained in consistency with the correlational] 



sin 2 e w (m z ) = 0.2324 - 10~ 7 x \[m t /GeV) 2 - 143 2 } ± 0.0003 (43) 

We have converted this correlation into a relation between sin 2 9w and tanfl, using the 
relation of m{ xd and mt , i.e. 

mt = m{ xd sin(3 (44) 

We first search for the tan[3 's satisfying the above correlation. Then, this range is 
further constrained by the requirement hb = h To at Mjj. In the present work, we have 
concentrated in the small tan(5 scenario, i.e. when ht 3> h\, T and we comment for the 
large tanfi case later. 

At the next stage, we introduce the Dirac neutrino RGE and run all of them together 
from Mjj down to the right handed neutrino scale Mat. We compare the predictions 
with those of the previous running (i.e. when there is no right-handed neutrino in 
the theory) and calculate the deviation from the equality of the r — b unification for 
the same inputs at Mrj. Below mjv we add the RGE for the effective light neutrino 



6 For example, we have found that domain walls may give structure at medium and large 
scales if, cither they are unstable, or the minima of the potentials of the relevant scalar field 
appear with different probabilities p0(]. 
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Majorana mass matrix. We assume that we are in a diagonal basis, so we can run the 
three eigenvalues of Mn independently. 

Let us start with the low tan/3 regime, assuming an effective SUSY scale M% Sj < 
ITeV. We vary ajj in a range close to a central value ^ and Mjj close to 10 16 GeV. 
Our first observation is that the introduction of an intermediate scale where the right 
handed neutrino gets a mass, shifts slightly the range of the parameter space for which 
unification is possible. For example, assuming ht G ~ 3, i.e., close to its infrared fixed 
point, and assuming a unification point ranging in Mq ~ (1.2 — 2.2) x 10 16 GeU, with 
k, (23.81 - 25.64), the effective scale ranges from M e s ff = (lOO)Gey to ITeV. 
Introducing the right-handed neutrino, we find M e s ff > UOGeV. Some particular 
cases with the corresponding low energy predictions are shown in Table 1. 



M N 
GeV 


1 

an 


M„ 

10 lb GeV 


M e / f 


l 


sinl 


«3 


10 16 


23.81 


2.18 


100 


127.9 


0.2325 


0.121 




23.81 


2.41 


110 


128.96 


0.2320 


0.123 




24.39 


1.97 


221 


128.09 


0.2320 


0.120 




25.00 


1.46 


493 


127.98 


0.2321 


0.118 




25.64 


1.08 


1212 


127.82 


0.2319 


0.116 


10 11 


23.81 
23.81 
24.39 
25.00 
25.64 


2.18 
2.41 
1.97 
1.46 
1.08 


110 
122 
270 
493 
1212 


127.83 
127.90 
127.89 
128.05 
127.89 


0.2323 
0.2318 
0.2315 
0.2321 
0.2320 


0.122 
0.124 
0.122 
0.118 
0.115 



We should point out that, the presence of vr in the spectrum has little effect in the low 
energy values of a em , sin 2 6w , «3 parameters. Moreover, for the above initial conditions 
the sin 2 6w — mt correlation, restricts tanj3 very close to unity tan/3 < 2. Of course, 
the biggest effects from the vr threshold are found in the b — r unification. 

For values in the perturbative regime of the top Yukawa coupling, ht G , at Mqut we 
find it impossible to obtain the correct mi,,m T masses starting with hb = h T at Mgut, 
even if the neutrino threshold is as high as Mjy = O(10 15 GeV). For example, using 
ht G ~ 3.2, (i.e. very close to its non-perturbative regime) and hb w h TQ w .0125, one 
can hardly obtain a running mass for the bottom m,b(mb) ~ 4.5GeV while the upper 
experimental bound is 77^(771,5) < AAGeV . However, the solution of the solar puzzle 
needs M N < 10 13 G e y. 

Therefore, in the following we do a complete two loop analysis and calculate the exact 
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deviation from the b — r - universality for a reasonable range of the scale Mjy. In our 
approach we first require the r - mass to be 1.749 ± 0.001 at mz- Then we search 
for the correct bottom mass and top mass as well as the required tan(5. We choose 
the biggest possible coefficient for which we have a solution, which corresponds to a 
bottom mass AAGeV . The variation of this coefficient as a function of M/v is plotted 
in Fig 1 (h T0 = 0.012), for h ta = 3.2 and 2.0 denoted in the plot by stars and crosses 
respectively. For the rest of the input parameters we take: 

M e J f = 5UGeV, Mjj = 1.46 x 10 16 G e y, a v = ^- (45) 

We see, in agreement with the qualitative analysis, that for this parameter range and 
small ht G it is not possible to obtain solutions for the b — r ratio at unification being 
unity. The larger the Yukawa coupling for the top, the lower the neutrino scale for 
which we find solutions. 

It is useful to compare the mass and other parameter predictions with respect to those 
obtained without the inclusion of vr. As has been pointed out in the previous section, 
in the presence of vr correct predictions for b — r masses can be restored either by 
increasing h TQ or by shifting hf, to smaller values at Mjj with a simultaneous change 
in the (3 - angle. In this latter case, we show in Fig 2 the curves corresponding to 
the values of tan(3 as a function of Mjv, needed to compensate for the decrease of the 
bottom mass. We find that as Mat decreases there is a large effect on tan(3, which 
drops for the two different choices for the top Yukawa coupling, from a common value 
of 1.35 at M N ~ 10 16 GeV to 1.02 and 1.13 for h ta = 3.2 and 2.0 respectively, at 
Mjv ~ 10 11 GeV. This, combined with the running of the top Yukawa coupling to the 
fixed point (Fig 3), implies that we expect in this case a decrease in the top mass, as 
the qualitative description of the previous section has indicated (Fig 4)- The larger 
the initial value of the top Yukawa coupling and the smaller the initial value of h T , the 
biggest the effect on tan(5 through the running, and the larger the effect on the top 
mass. 

In Fig. 5 we see the effect of Mjv on l/a em , which increases slightly as Mjy drops. At 
the same time, sin 2 9 w also increases slowly (Fig 6), while the strong coupling decreases 
(Fig 7). In all cases the effect is much smaller than the errors on these quantities, 
however it is enough to eliminate some of the solutions that were in the border of the 
allowed region. 

We would like to stress that, in the case where the hb is the same while h To is slightly 
increased to compensate for the reduction caused by there is no need to change the 
angle [3. For this reason there is no significant effect on the top mass, which preserves 
its value as in the standard case. 

Finally, in Fig 8 we plot the light r-Majorana neutrino mass versus h ta coupling, for 
three different values of the heavy Majorana scale Mjv = 10 12 , 10 13 and 10 14 GeV. 
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This analysis can be applied also in the case of the large tan(3 regime. However in this 
case there are important corrections to the bottom mass from one-loop graphs involving 
susy scalar masses and the \i parameter. These graphs might induce corrections to 
even of the order of (30 - 50)%. 

In view of the considerable uncertainties involved we have not extended the numerical 
analysis to cover this case. 



7 Conclusions 



In this paper we have discussed the implications for low energy physics of right-handed 
neutrinos with masses at an intermediate scale Mjq. For Mjy ~ 10 12 ~ 13 GeV (required 
to give a r neutrino mass 0(leV) to serve as a hot dark matter component) a 10% 
deviation from b — r mass equality at the GUT scale is needed to give an acceptable 
value for the ratio mj, / m T as measured in the laboratory. We showed that it is possible 
to retain the m® = m® GUT prediction of the (3, 3) - elements of the corresponding 
mass matrices provided there is sufficient mixing in the charged lepton mass matrix 
between the two heavier generations. The scenario we propose can be realised in a 
simple extension of the standard symmetry of electroweak interactions to include a 
U(l) family symmetry. Consideration of the implications of this symmetry for neutrino 
masses shows that the large mixing implied allows for a simultaneous explanation of 
the atmospheric neutrino problem and the solar neutrino problem. This complements 
our previous discussion of solutions to the solar neutrino deficit while having a neutrino 
mass in the range needed to fit the COBE measurements in a hot plus cold dark matter 
universe. We have also presented detailed numerical solutions of the renormalisation 
group equations for the case of a heavy right-handed neutrino to support the analytical 
analysis. 
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Fig. 1 



Tanb versus neutrino scale 
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Running top Yukawa coupling 
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Fig. 3 



Top mass for decreased bottom coupling 
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Fig. 4 
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127.99 



Inverse electromagnetic coupling versus neutrino scale 



127.97 - 



l/a t 



127.95 



127.93 




26 28 30 32 34 36 

Log(M N ) 



Fig. 5 



sin 6 W versus neutrino scale 
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Fig. 6 
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Strong coupling versus neutrino scale 
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Fig. 7 



Light neutrino mass (in GeV) 
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